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The filamentary structures and associated electromagnetic fields of a relativistic electron beam 
have been studied by three dimensional particle-in-cell (PIC) simulations in the context of Fast 
Ignition fusion. The simulations explicitly include collisions in return plasma current and distinctly 
examine the effects of beam temperature and collisions on the growth of filamentary structures 
generated. 
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The transport of laser generated electron beams in 
overdense plasmas is an important topic to study spe- 
cially in the context of Fast Ignition (FI)i. PIC simula- 
tions and calculations have shown that the typical energy 
of electrons in these beams is of the order of a few MeV, 
which implies Giga- Ampere currents needed to transport 
Peta-Watt energy fluxea^ii^. However, the transporta- 
tion of beam currents greater than the Alfven current 
limit is not possible as the self-generated magnetic field 
associated with these currents bends electron trajecto- 
ries strongly thereby preventing the forward propagation 
of the beam. In vacuum, the Alfven current limit is 
laif = 17 j(3k A, where 7 = 1/^1-/3^ is the relativis- 
tic Lorentz factor associated with the beam. A current 
greater than this limit can only propagate as a charge 
neutralized current. This sort of situation occurs in the 
FI scenario, where the charge of the relativistic electron 
beam is compensated by the return plasma current and 
the propagation of the electron beam is possible. But this 
configuration, where two opposite streams of plasma cur- 
rents are counterpropagating, is susceptible to the elec- 
trostatic two-stream instability and the transverse fila- 
mentation instability, often referred as the Weibel insta- 
bility. Due to these instabilities, the electron beam gets 
splitted into many filaments^i^i^. This filamentation can 
significantly affect the beam energy deposition at the cen- 
ter of the fuel target and thus is one of the central issues 
for the FI concept. Recently, Honrubia et alX have car- 
ried out 3D hybrid simulations of the fast electron trans- 
port and resistive beam filamentation in inertial fusion 
plasma. Califano et al^ have reported three dimensional 
magnetic structures generated due to the Weibel instabil- 
ity of a relativistic electron beam in coUisionless plasmas. 
Three dimensional PIC simulations of the Weibel insta- 
bility in astrophysical context have also been reported^. 
Moreover, the evidence of Weibel- like dynamics and the 
resultant filamentation of electron beams have been ob- 
served experimentallyi^. 

The role of collisions in the Weibel instability is rather 
non-trivial and has been studied analytically in linear and 



quasilinear regime of the instabilityi^. The effect of col- 
lisions on the nonlinear regime of the instability is rather 
unaccessible to analytical calculations and can best be 
studied by PIC simulations. However, till now the full 
scale PIC simulations on the filamentation of the elec- 
tron beams have mostly concentrated on the coUisionless 
Weibel instability and have not explicitly included the 
effect of the collisions in the return plasma current. At 
the same time, collisions in the return plasma current 
play an important role on the Weibel instability. The 
eigenmodes of a collisional plasma have been computed 
recently, where it was shown that frictional nature of col- 
lisions can change the dynamics of the eigenmodes of the 
system significantlj*i^. In the presence of collisions the 
Weibel instability can not be suppressed by the trans- 
verse beam temperature as proposed initially^^-. A more 
focused discussion about the collisional effects on the 
Weibel instability has also been addressed recently where 
it was shown that even small collisions in the background 
plasma can trigger the Weibel instability of a warm elec- 
tron beam that would be otherwise stable in coUisionless 
plasmaii^. 

Motivated by these findings, we carry out three- 
dimensional PIC simulations on the transport of a rel- 
ativistic electron beam in the FI plasmas by explicitly 
including the effect of collisions in the return plasma 
current. The rationale behind choosing the 3D geome- 
try is twofold: the first reason is to study the realistic 
geometry and the second is to show that in this geome- 
try, the Weibel instability can't be suppressed due to the 
transverse beam temperature. It has been shown before 
analytically that in the 3D geometry the Weibel instabil- 
ity survives even with high transverse beam temperature 
due to the presence of an "oblique mode" , which is the 
manifestation of coupling between the Weibel instabil- 
ity and the two-stream instability^. We conjecture the 
existence of the "oblique mode" due to the anomalous 
coUisionality of the background plasma. It arises due 
to the turbulence of electrostatic fields generated by the 
two-stream instability. We choose four cases for the sim- 
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ulations (A) cold electron beam in a collisionless back- 
ground plasma, (B) cold electron beam in a coUisional 
background plasma, (C) warm electron beam in a cold 
collisionless background plasma, and (D) warm electron 
beam in a coUisional background plasma, thus highlight- 
ing the influence of these physical processes separately. 
Hence, for the first time we have studied the filamenta- 
tion of the relativistic beam investigating the role of each 
physical process clearly in three dimensional particle-in- 
cell simulations. 

In our simulations, the electron beam propagates along 
the negative .Z— axis with the relativistic velocity V(b,z) ^ 
^(fe.a;) J ''^(b,y) ■ The bulk cold background plasma is repre- 
sented by ambient plasma electrons, while the plasma 
ions are considered as a fixed charge-neutralizing back- 
ground with the density uq = ni, + Up, where Ub and Up 
are the densities of the electron beam and the background 
plasma, respectively. Initially the beam current is fully 
compensated by the return plasma current. The beam 
density is much smaller than the background plasma elec- 
tron density, i.e. rih ^ Up, a, usual scenario of the FI. The 
spatial dimensions of the simulation domain L is large in 
comparison with the electron skin depth i.e. L ^ Xs, 
where As = c/ujpe; c and tOpe are the velocity of light 
in vacuum and the total electron plasma frequency, re- 
spectively. The quasi-neutrality is maintained over all 
simulation time as the field evolutions due to the Weibel 
instability occur on a time scale slower than the plasma 
electron frequency At >> l/wpe- The coUisional pro- 
cesses are simulated with a newly implemented collision 
module in the relativistic PIC code Virtual Laser Plasma 
Laboratory (VLPL)i^. 

The 3D simulation box has spatial dimensions X xY x 
Z = {20Xs X 20As X 20As). The 3D simulation domain is 
sampled with the mesh of 160 x 80 x 20 cells. All simula- 
tions are performed with 64 particles per cell and with a 
grid size much smaller than the background plasma skin 
depth (hx, hy,hz << Xs). The density ratio between the 
beam and plasma electrons is Up/ni, — 9, whereas the 
beam and the background plasma electrons have veloc- 
ities Vb = 0.9 c and Vp = 0.1c. The evolution of the 
field energies for every component Fi of the fields E and 
B is recorded at every diagnostic step of the simulation 
summed over all the grid cells as 

/ dxdydz = {eFiirieCUjpef Khyhz, (1) 
Jv 

where, h^, hy and hz are the spatial grid sizes along x, 
y and z and eFi/meCLOpe represents the relativistic field 
normalizations. The relativistic electron beam has a tem- 
perature of Tf, « 70 keV and the ambient plasma collision 
frequency is assumed to be Vei/ojpe — 0.15 for these sim- 
ulations. The background plasma is always cold initially 
while the beam electrons do not face any collisions. 

Fig. [T] shows the structure of the beam filaments dur- 
ing the nonlinear regime at time T — 16(27r/a;pe) for 
the case (A). The filaments transverse spatial sizes are 








FIG. 1: (color online) 3D structure of the beam filaments 
generated in the PIC simulation at a time T = l&{2n / uipe) 
in the nonlinear regime for the case (A), where the electron 
beam is cold and the plasma is collisionless. The filaments 
are tiny and their transverse dimension is comparable to the 
skin length of the cold background plasma. 



comparable to the plasma skin depth As while their lon- 
gitudinal spatial extent is larger than the plasma skin 
depth. This could be understood as follows: in the linear 
regime of the beam filamentation, the transverse extent 
of the filaments is governed by the Weibel instability and 
its growth rate maximizes around the plasma skin depth 
in /c-space. However, the longitudinal extent is deter- 
mined by the two-stream instability which maximizes at 
wavelengths much larger than the plasma skin depth. As 
a result, the filaments have small spatial extent in the 
transverse direction and elongated in the propagation di- 
rection. Later on, during the nonlinear stage of the insta- 
bility, the filaments merge together and their radii grow. 
The growth rate of the filamentation instability depends 
on the beam-to-plasma densities ratio, and is expressed 
as F = (■\/3/2'*/^) (nb/np7h)^^'^, where 7b is the relativis- 
tic Lorentz factor of the electron beaniii. The Weibel in- 
stability generates very strong quasi-static magnetic field, 
which can be seen from the surface plot of the transverse 
(Bj;, By) and longitudinal {Bz) magnetic field B^. shown 
in Fig. [21 Evidently, the amplitude of the transverse com- 
ponents is larger at least by one order of magnitude com- 
pared to the longitudinal component. The contour lines 
on the bottom surface in this figure show that each of 
these tiny filaments is surrounded by strong axial mag- 
netic fields. We have plotted the temporal growth of such 
field energies for all four simulation cases later in Fig. [B] 

The simulation case (B) is shown in Fig.[3l We observe 
that the transverse spatial width of the filaments is larger 
while the longitudinal extent remains unchanged. Due 
to collisions, the growth rate of the Weibel instability 
saturates at large plasma skin depth in fc-space. At later 
times, the filaments merge together and their transverse 
spatial extent grows even larger, which is evident from 
the simulation results. 
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FIG. 2: (color online) Illustration of the transverse {Bx, By) 
and longitudinal (Bz) magnetic fields at the same time as 
Fig. [T] The contour lines on the bottom surface gives a 2D 
projection. B fields are normalized as explained in the Eq. ([T]). 
Magnitude of the transverse components are larger almost by 
an order than the longitudinal component. 




FIG. 4: (color online) Beam filamentation persists in the case 
(C), where the electron beam is warm with Tb « 70keV prop- 
agating in a collisionless plasma, at time T — 10{2-K/ujpe) 
in the nonlinear regime. The evacuation of the background 
plasma is not complete due the strong thermal effects. 








FIG. 3: (color online) 3D structure of the beam filaments 
from the PIC simulation in a time T = 8{2-K/ijpe) for case (B), 
cold electron beam in collisional plasma. The filament sizes 
are larger than collisionless case due to collisional dissipation. 



The filament structures of the simulation case (C), 
where the electron beam has a transverse beam tem- 
perature Tb w 70 keV, is shown in Fig. 0] at time T — 
10{2TT/ujpe). The Weibel instability could be suppressed 
by the high transverse beam temperature. However, in 
3D geometry the coupling of the Weibel and the two- 
stream instabilities can manifests itself into the so called 
oblique mode, which persists even when the beam tem- 
perature is high. We may also consider the occurrence 
of the beam filamentation in this case due to "anoma- 
lous plasma coUisionality" as shown in Refji^. The two- 
stream mode actually generates electrostatic turbulence 
in the plasma. Stochastic fields associated with this tur- 



FIG. 5: Structure of the electron beam filaments derived from 
the PIC simulations for the case (D), where the warm electron 
beam is propagating in the collisional background plasma. 
The filaments are larger and total evacuation of the back- 
ground plasma is prominent. 



bulence scatter the beam and plasma electrons and lead 
to an effective coUisionality in the return current. This 
anomalous effect revives back the Weibel instability. The 
phenomena of anomalous resistivity of plasma generated 
due to the turbulence of electric and magnetic fields have 
been noted earlier also^.. The allowance for finite back- 
ground plasma temperature decreases the growth rate 
of the instability. The interplay of collisions in different 
regimes of beam and background plasma temperatures on 
the Weibel instability can be found in Ref.^^. Moreover, 
it maybe worthwhile to mention here that this filamen- 
tation is lacking complete evacuation of the background 
plasma electrons contrary to the previous cases, presum- 
ably due to the strong thermal effects. 

Fig. O shows the beam filaments of the last simulation 
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FIG. 6: (color online) Time evolution of Weibel field ener- 
gies {E^ and B^) for the different simulation cases: (A) cold 
electron beam in a coUisionless background plasma, (B) cold 
electron beam in a coUisional background plasma, (C) hot 
electron beam in a coUisionless plasma and (D) hot electron 
beam in coUisional background plasma. The time is normal- 
ized in units of 27r/t<;pe and the E and B fields are normalized 
as (eE/meOjpec) and (eB/meWpec), respectively. 

case (D), where the electron beam is hot and the back- 
ground plasma is coUisional. Here again we see strong fil- 
amentation and the spatial width of the beam filaments 
are also larger due to collective coUisional plasma effects. 
Moreover, we have observed almost complete evacuation 
of the background plasma in this case. 

To get further insight into the growth of the fields asso- 
ciated with the filamentation instability, we have studied 
the the evolution of electric and magnetic field energies, 
shown in Fig. \6\ for the four different simulation cases. 



The vertical axes in Fig. [6] represents the normalized field 
energies in logarithmic scales whereas the horizontal axes 
are for time scaled in 27r/ujpe. In each of these cases, 
we see a stage of linear instability, where the field ener- 
gies build up exponentially in time and it is followed by 
the nonlinear saturation of the instability. Afterwords, 
in the nonlinear stage of the instability, the filaments 
merge rapidly with each other due to the magnetic at- 
traction and the field energies saturate. One may note 
that the simulation cases (B), (C), and (D) display simi- 
lar trend of field energy build up. This similarity is very 
interesting for the case (C), where the electron beam is 
warm and has the sufficiently high temperature to kill 
the Weibel instability. Yet, we still see the filamentation. 
Clearly it happens due to the presence of the two-stream 
or "oblique mode" . Based on the similarity between cases 
(B) and (C), we attribute the filamentation of the beam 
to the 'effective plasma coUisionality' of the beam system 
as discussed before in the manuscript as welU^. 

In summary, we have studied the filamentation in- 
stability of a relativistic electron beam in Fast Ignition 
plasma by three dimensional particle-in-cell simulations, 
distinctly examining the effects of high beam tempera- 
ture and plasma collisions. The important result of this 
study is that the beam temperature does not suppress the 
filamentation instability even in absence of collisions in 
beam plasma system. An explanation on the persistence 
of the Weibel instability in 3D geometry is offered. It is 
attributed to the anomalous coUisionality of the beam- 
plasma system due to the two-stream mode. The non- 
linear regime of the filaments merging and the growth of 
field energies due to plasma collisions have been studied. 

This work was supported by the DFG through TR-18 
project and the US Department of Energy grants DE- 
FG02-04ER41321 and DE-FG02-04ER54763. 
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